How magic is the magic ^^Ni nucleus? 
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We calculate the B{E2) strength in ^*Ni and other nickel isotopes using several theoretical ap- 
proaches. We hnd that in ®*Ni the gamma transition to the Hrst 2^ state exhausts only a fraction 
of the total B{E2) strength, which is mainly collected in excited states around 5 MeV. This effect is 
sensitive to the energy splitting between the fp shell and the orbital. We argue that the small 
experimental B{E2) value is not strong evidence for the double-magic character of ®*Ni. 



The appearance of shell gaps associated with magic 
nucleon numbers is one of the cornerstones of nuclear 
structure. The presence of magic gaps allows one, for 
example, to determine the single-particle energies and 
the residual interaction among valence nucleons, provid- 
ing essential input for nuclear models. Magic gaps offer 
a natural way of performing truncations in microscopic 
many-body calculations. Magic nuclei also play an es- 
sential role in the two major nucleosynthesis networks 
(s- and r-process) that produce the majority of nuclides 
heavier than mass number A ^ 60. 

The doubly-magic character of ^^Ni (Z=28, iV=40) 
was suggested in the early eighties [1,2] and tested ex- 
perimentally [4-7]. The proton number Z—28 in the 
nickel isotopes is magic. In the neutrons, the sizeable en- 
ergy gap at A^=40 separates the pf spherical shell from 
the gg/2 intruder orbit. However, this spherical subshell 
closure is not sufficiently large to stabilize the spheri- 
cal shape. [Experimentally [3], *°Zr (N^Z^AO) behaves 
like a well-deformed rotor.] The current experimental 
evidence about the double-magicity of ^®Ni is controver- 
sial [8]. On the one hand, ^*Ni does not show a pro- 
nounced irregularity in the two-neutron separation ener- 
gies, as is expected for a magic nucleus. On the other 
hand, the lowered position of the Oj level, the slightly 
elevated energy of the first 2+ state, and the quite small 
i3(i?2,0+g 2f) value are often interpreted as indi- 
cations for magicity. In fact, the B{E2) rate in ^^Ni 
(280 ± 60 e^fm^ [7]) is significantly smaller than that in 
the well-established double-magic nucleus '^^Ni (620±120 
e^fm^ [9]). 

As discussed in Ref. [10], the size of the N—4Q gap 
strongly depends on the effective interaction used, and 
it dramatically influences the quadrupole collectivity of 
the 7V=40 nuclei. While there is much discussion in the 
literature about the weakening of shell effects in neutron- 
rich nuclei (e.g., the magic gap 7V=28 seems to be eroded 
in drip-line systems; see Ref. [10] and references quoted 



therein), ^^Ni lies very far from the expected neutron 
drip line (expected to be around ®^Ni [11]) and one 
should probably not invoke "exotic" explanations when 
discussing the stucture of this neutron rich nucleus. 

It is the aim of this Letter to draw attention to the 
total B{E2;0+^ -> 2^ distribution in '^^Ni, which can 
hold the key to the question whether this nucleus is magic 
or not. We will argue that the transition to the first ex- 
cited 2+ state constitutes only a small part of the total 
B{E2) strength and that the total B{E2) strength de- 
pends sensitively on the size of the iV = 40 shell gap. 

To understand the structural difference between ^^Ni 
and ^^Ni (where the transition to the first 2+ state ex- 
hausts most of the total low-energy strength), we be- 
gin from qualitative arguments based on a simple inde- 
pendent particle model (IPM). Proton configurations in 
both nuclei are identical and correspond to a closed (77/2) 
shell. The neutron configuration in ^^Ni is the same as 
the proton one, (77/2)^, while in ^^Ni it corresponds to 
the closed {pf) shell. The first excited 2+ state can be 
viewed as a superposition of particle-hole (ph) excita- 
tions. In ^^Ni, the lowest proton and neutron ph exci- 
tations are identical (due to isospin symmetry) and they 
are of (Ip-lh) character. The 2'1 state can be roughly ex- 
pressed as l/\/2(<I>p(lp- Ih) + <I>„(lp- Ih)). Due to the 
parity change between {pf) and 39/2 orbits, the 2^ state 
in ^^Ni cannot have a Ip-lh neutron component. It can 
thus be written as a$p(lp — Ih) -\- (3^n{'^P — 2ft.), where 
|q;|^ + |/5|^=1- The weight of the neutron compo- 
nent decreases with the size of the iV=40 gap. As B{E2) 
transitions reflect a proton component in the wave func- 
tion, an increasing value of j3 reduces the B{E2) rate. 
Consequently, a smaller B{E2,0+^ 2^) value in '^'^Ni 
than in ^^Ni suggests that [al is noticeably smaller than 
l/\/2, implying that it is more favorable in ^^Ni to excite 
the pair of neutrons into the gg /2 orbital than to excite a 
single proton across the magic = 28 gap. As the pro- 
ton configurations in the ^^Ni and ^^Ni ground states are 
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the same, the IPM also suggests that a noticeable part of 
the (Ip-lh) proton excitations in ^^Ni should be mixed 
into excited states. Of course, the neutron amplitude 
hence the B{E2) rate, should also depend on the residual 
interaction. Therefore, to check whether the simple IPM 
argument also holds in the presence of a realistic resid- 
ual interaction, we performed calculations in three dif- 
ferent theoretical models: the Shell Model Monte Carlo 
(SMMC), the Quasi-particle Random-Phase Approxim- 
mation (QRPA), and a large-scale diagonalization shell 
model (SM). 

The SMMC approach allows the calculation of nuclear 
properties as thermal averages, employing the Hubbard- 
Stratonovich transformation to rewrite the two-body 
parts of the residual interaction by integrals over fluctu- 
ating auxiliary fields [12]. We performed SMMC studies 
of the even-even nickel isotopes between ^^Ni and ^^Ni in 
the complete {fp){gds) configuration space for both pro- 
tons and neutrons. The single-particle energies were de- 
rived from a Woods-Saxon potential appropriate for ^^Ni, 
placing the important levels at excitation energies (in 
MeV) of 4.3 (P3/2), 6.4 (f,/^), 6.6 (pi/2), 9.0 (39/2), 13.0 
(^5/2) relative to the /7/2 orbital. We employed the same 
residual interaction of the type pairing-plus-quadrupole 
as in Ref. [15], which allowed us to avoid the sign prob- 
lem in SMMC calculations [13]. We checked that our in- 
teraction gives a reasonable description of the collective 
spectrum of ^^Ge and ^^Ni, and that center-of-mass con- 
taminations are small and do not affect our results for 
quadrupole excitations. The SMMC calculations were 
performed at temperature T=0.33 MeV (corresponding 
to 96 A/3 'time slices'), which, for even-even nuclei, is 
sufficiently low to cool the nucleus to the ground state. 
We checked this for ^*Ni and found variations of the var- 
ious quadrupole expectation values of less than 3% by 
slightly increasing the temperature to T=0.4 MeV. The 
Monte Carlo integrations used between 1000 and 4000 
samples. 

In the SMMC approach, the total B{E2) strength is 
approximated by the expectation value B{E2) oc (Q^), 
where the quadrupole operator is defined by Q = epQp + 
SnQn, with (5p(„) = Y,i rfY2{9i,4>i); the sum runs over aU 
valence protons (neutrons). The effective charges ep,e„ 
account for coupling to the states outside our model 
space. We adopt in the following the standard values 
Bp = 1.5, e„ = 0.5 [14]. For the single-particle wave 
functions, we adopt the harmonic oscillator states with 
the oscillator length b = l.OlA^/^ fm. At low temper- 
atures, the total B{E2) strength obtained in SMMC 
approximates the expectation value (O^g [Q^jO+g ) — 

E/l(0+s.lQ|2;)|' « EfB{E2;0+,, 2+). Therefore, 
the total SMMC strength corresponds to the summed 
B{E2) strength to the excited 2+ states within the as- 
sumed configuration space. 

The total SMMC B{E2) strengths are plotted in Fig. 1. 
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FIG. 1. Comparison of the total SMMC B{E2) values for 
even-even nickel isotopes (solid circles) with the experimental 
B(£;2,0+s. 2+) rates (open circles, from Ref. [16]). The 
solid squares represent the SMMC B{E2) values obtained in 
the SMMC variant of calculations in which all the gds sin- 
gle-particle energies are shifted up by 1 MeV. 

closely up to ^^Ni. For these nuclei, it is well known 
from electron scattering experiments that most of the 
B{E2) strength resides in the transition to the first 2+ 
state (see, e.g., Ref. [17]). For the isotopes approach- 
ing TV = 40, SMMC predicts a significantly larger to- 
tal B{E2) strength than observed in the first transition. 
For "^^Ni, the calculated total strength (-- 1250e2fm'*) is 
about five times greater than the measured transition to 
the first 2+ state. This is consistent with the fact that 
the centroid of the SMMC B{E2) strength, calculated 
from the respective response function, lies at ~ 5 MeV. 
This value is significantly higher than the energy of the 
first 2+ state and indicates that most of the calculated 
SMMC B{E2) strength in resides in excited states. 

In Ref. [7] it was pointed out that the neutron gg/2 or- 
bital plays a major role at = 40, thanks to cross-shell 
neutron pairing excitations. We confirm this finding. In 
our SMMC study we find an average occupation num- 
ber (n)=2.2 for neutrons in the 179/2 orbital (and 0.22 
in the c?5 /2 orbital which couples strongly to 179 /2 by the 
quadrupole force). These numbers are significantly re- 
duced (to 0.9 and 0.08, respectively) if one artificially 
shifts upwards all levels of the gds shell by 1 MeV, mak- 
ing the A^=:40 gap larger (see the SMMC variant of cal- 
culations in Fig. 1). As a consequence, the quadrupole 
moment of the neutron configuration gets reduced and 
the summed B{E2) value decreases. Although in this 
modified calculation the summed B{E2) value is smaller 
for ^*Ni than in the neighboring nuclei, our calculation 
still predicts most of the strength in excited states. We 
note that for both sets of single-particle energies, ^*Ni 
is characterized as magic, having a reduced B{E2) value 
compared to the neighbors. It is also worth mention- 
ing that the variations of along the isotope chain are 
rather small. The smallest value is obtained for ^^Ni, the 
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largest for ^^Ni; however, the variation is less than 9%. 
This shows again that proton (/7/2 P3/2,Pi/2, /5/2) ex- 
citations cannot be neglected in the reproduction of the 
E2 strength, as already noted in [7] , but the dominating 
variations in the B{E2) strength come from the neutrons. 

While the SMMC approach allows for the calculation 
of the summed strength in large model spaces, it is not ca- 
pable of making detailed spectroscopic predictions. For 
this reason, we have also performed QRPA and diago- 
nalization shell-model calculations. Our QRPA calcula- 
tions closely follow the formalism described recently in 
Ref. [18]. As a residual two-body interaction, we use the 
sum of an isoscalar and an isovector quadrupole force, 
and a quadrupole pairing force. For the single-particle 
levels below the N {Z) = 82 shell gap. we took those of 
the Woods-Saxon potential [20], and the unbound states 
were approximated by the Nilsson levels [21]. Guided by 
the experimental data (cf. Ref. [19]), the energy of the 
2p3/2 neutron state was shifted up by 1 MeV. For the 
strength of the isoscalar quadrupole force, we adopted 
the self-consistent value multiplied by 0.8, and for the 
isovector force we took xt=i = ~^^tt \^^\- The renor- 
malization factors of the pairing gaps are 0.8 (neutron) 
and 0.9 (proton). The bare charges are used for calcula- 
tion of B{E2; — > 2+), consistent with the large model 
space. The results are shown in Fig. 2. 
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(a) Comparison of the QRPA energies 
summed B{E2) strength (filled squares) and 
J. 2^) values with the experimental data. The 
summed B{E2) strength includes all the transitions up to an 
excitation energy of 9 MeV. 



FIG. 2. 
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Our calculations nicely reproduce the observed trend of 
the Er^+ energies in Ni isotopes, including the pronounced 

rise at ^^Ni and ^*Ni. The QRPA calculations also give 
a reasonable description of the B{E2, 0+g ^2^) values, 
with the maximum around ^^Ni and the strong decrease 
towards ^*Ni. The structure of the lowest 2+ QRPA 
phonon in ^^Ni is dominated by neutrons (90%). Im- 
portantly, our QRPA calculations confirm that most of 
the B{E2) strength in ^^Ni resides in excited states, in 
contrast to ^^Ni, where the B{E2) strength is exhausted 
by the transition to the first 2+ state. These arguments 
are demonstrated again in Fig. 2, which displays the 
summed B{E2) strength (filled squares) and in Fig. 3, 
which shows the predicted B{E2) strength function. It 
is seen that around neutron number iV = 40, a signifi- 
cant quadrupole strength develops at excitation energies 
around 4 MeV. 



CU 0.04 



(N 

T 

+ 

o 
w 



Da 



°Ni - 



"Ni - 



°Ni - 



°Ni 



°Ni 



"Ni 



J_ 



E (MeV) 

FIG. 3. Distribution of the B{E2;Q+^, 2+) strength for 
even-even Ni isotopes calculated in the QRPA method. 



Finally, our arguments have been tested and confirmed 
in large-scale diagonalization shell-model calculations of 
the B{E2) strength distribution in 66,68]vgj^ adopt the 
same valence space and the effective interaction fpg cm- 
ployed in Ref. [7]. This valence space consists of a ^*Ca 
core (more precisely, a '*''Ca core with eight f-j/2 frozen 
neutrons), the f'j/2, ^3/21 P1/21 ^-nd active orbitals 
for protons and the P3/2, P1/2, and c/9/2 active or- 

bitals for neutrons. The SM calculations describe well 
the behavior of the 2+ energies and B{E2) rates for the 
Ni isotopes ranging from N=28 to N=40 [7]. The B{E2) 
strength distribution shown in Fig. 4 has been calculated 
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2"*") strength in 
®®'^*Ni calculated in the diagonalization shell model. 



the B{E2) strength should reside in excited states. Fur- 
ther experimental investigations, including the g-factor 
measurement in ^^Ni, are certainly called for. 
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at a truncation level which considered up to a total of 7 
particle excitations from the /y /2 orbital to the upper fp 
shell for protons and from the upper fp shell to the gg/2 
orbital for neutrons. For ^^Ni, the calculation predicts 
a B{E2) value of about 280 e^fm^, which nicely agrees 
with the experimental value. However, the transition to 
the first 2+ state exhausts only the smaller fraction of the 
total shell-model B{E2) strength, which we calculate as 
900 e^fm^. Most of the calculated strength resides at ex- 
citation energies around 5-6 MeV. The wave function of 
the first 2f state involves mainly (2p-2h) neutron exci- 
tations; the proton configuration is a mixture of Op-Oh 
(50%), Ip-lh ((25%), and 2p-2h (20%) components. The 
SM ^^Ni ground state corresponds to a closed-shell con- 
figuation plus a 35% admixture of 2p-2h neutron excita- 
tions. The results for ^^Ni nicely confirm the QRPA pre- 
diction: the transition strength is gradually shifted from 
the first excited state to higher energies as approaching 
N=AO. 

In summary, we have performed microscopic calcula- 
tions of the i?(i?2;0+g 2+) strength distribution in 
^*Ni, and in other even-even nickel isotopes. Our main 
finding is that a significant portion of the B{E2) strength 
in ^^Ni resides in excited states around 5 MeV, and that 
the small observed B{E2, O+g 2^) value is not nec- 
essarily an argument for a shell closure at iV = 40. In 
fact, we argue that this transition rate is quite sensitive 
to the energy splitting between fp shell and 59/2 orbital, 
and that its smallness might indeed be an indication for a 
rather small gap. In short, the small B{E2) rate refiects 
the fact that the lowest 2+ state in ^^Ni is primarily a 
neutron excitation (cf. Ref. [18] for a similar discussion 
for ^^^Te). If our arguments are correct, then most of 
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